THE SUPPORT OF THE LIMIT DISTRIBUTION OF OPTIMAL 
RIESZ ENERGY POINTS ON SETS OF REVOLUTION IN 



J. S. BRAUCHART*, D. P. HARDINf, AND E. B. SAFFf 



Abstract. Let A be a compact set in the right-half plane and r(A) the set in 
obtained by rotating A about the vertical axis. We investigate the support 
of the limit distribution of minimal energy point charges on r(j4) that interact 
according to the Riesz potential l/r", < s < 1, where r is the Euclidean 
distance between points. Potential theory yields that this limit distribution 
coincides with the equilibrium measure on T{A) which is supported on the 
outer boundary of T{A). We show that there are sets of revolution T{A) such 
that the support of the equilibrium measure on r{A) is not the complete outer 
boundary, in contrast to the Coulomb case s = 1. However, the support of the 
limit distribution on the set of revolution r{R + A) a,s R goes to infinity, is 
the full outer boundary for certain sets A, in contrast to the logarithmic case 
{s = 0). 



1. Introduction 

The discrete energy problem for Riesz kernels fcs(x):=|x|^*, s > 0, on compact 
sets K in M'^ is concerned with finding iV-point systems in K in the most-stable 
equilibrium; that is, that minimize the s-energy 

1 N N 

among all iV-point sets Xf^ :— {xi, . . . ,XAr} c K, where | • | denotes Euclidean 
distance. The existence of such configurations follows from both the lower semi- 
continuity of the Riesz kernel fc^, s > 0, and the compactness of K. Even in the 
case that K is the unit sphere in R'^, explicit examples of such point sets are known 
only for a few values of N . For approximate physical models of configurations of 
minimal energy points for large N on the sphere as well as toroidal surfaces, see 
[5, 6]. 

The A^-point system defines a discrete measure Ai(-'i^7v):=(l/-^) X^xex ^x, 
by placing the charge at every point x G Xj^. In this paper we investigate 
the support of the limit distribution (limit in the weak-star sense as A'^ — s- oo) of a 
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sequence of measures /x(X^), N >2, induced by minimal energy point configura- 
tions X'^ on sets of revolution r(^) in B? obtained by revolving a compact set A 
in the right-half plane about the vertical axis. 

If < s < dimr(^) (the Hausdorff dimension of r(A)), classical potential theory 
for the Riesz kernel kg (cf. [12]) can be used to study this problem. In this case, 
the limit distribution (as TV —t oo) of optimal A''-point configurations is given by 
the equilibrium measure ^1s,t{A) that uniquely minimizes the continuous energy 

^■sM — yy ^6(x-y)d/x(x)d/i(y) 

over the class M{V{A)) of (Radon) probability measures /i supported on r(^). 
(For example, when T{A) is the unit sphere in R'^ the equilibrium measure is 
the normalized surface area measure on 5^.) 

The probability measure ijls,t(A) is characterized by the following variational 
principle [12, Ch. II]: For T{A) there exists a constant Vs = Vs{T{Ay) such that 

(1.2) Us"'^^^^ > Vs "approximately everywhere" on T{A), 

(1.3) Us"'^^^^ < Vs everywhere on the support of Hs,r(A)- 
Here Us"'^^"^^ denotes the equilibrium potential 

[/^•^<-*'(x):= j fc,(x-y) dM«,r(A)(y), x G R3_ 

The constant Vs is the infimum of the energies of (Radon) probability measures 
supported on T{A), that is Vg — Is[/^s,r(A)]- The reciprocal of Vs is called the 
s-capacity of the set T{A), it is denoted by capgr(A). The term "approximately 
everywhere" means that the property holds everywhere with the possible exception 
of a set of s-capacity zero. It follows from (1.2) and (1.3) that Ug''^^^^ = Vs 
approximately everywhere on the support of /is,r(A)) which provides an integral 
equation for the equilibrium measure on its support. Knowing this support is 
therefore an important step in the determination of iis,r{A)- 

We remark that Fabrikant et al. [9] provide a method for finding the density p 
of a signed charge distribution for a prescribed /cg-potential distribution on certain 
surfaces of revolution in R^. However, their methods do not apply, for example, 
to the torus and, more importantly, the distribution they obtain need not be non- 
negative. For the analysis of charge distributions in the Coulomb case (s = 1) on 
circular or ellipsoidal "slender toroidal surfaces", see Cade [7] and Shall [15]. 

Several important properties of the Riesz equilibrium measure ^s,k for a com- 
pact set K of positive s-capacity are summarized in the previously cited book of 
Landkof. Adopting the same notation, we let Goo denote the unbounded connected 
component of the complement of K. The boundary S of Goo is called the outer 
boundary of K. Furthermore, let K be "the set of all points of K each neighbor- 
hood of which intersects K in a set of positive s-capacity" ([12, Ch. II, no. 13]). In 
the case 1 < s < dimii', the First Maximum Principle yields that supp ^s,k 3 S. 
In particular, if s = 1, then supp/i.5_;^ = S. For s < 1, it follows from the supcr- 
harmonicity of the kernel ks that the equilibrium measure is concentrated on the 
outer boundary S of K. In [11] Hardin, Saff, and Stahl proved a stronger result 
for the logarithmic case (limit as s 0"*"): For any compact set A in the interior 
of the right half-plane IHI+, the limit distribution of minimal energy point charges 
on T{A) that interact through a logarithmic potential log(l/|x — y|) is supported 
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on its "outer-most" portion only. The "outer-most" part of a torus, for example, is 
the set of revolution generated by rotating the right semi-circle about the vertical 
axis. Numerical experiments (cf. [17] and Section 6) suggest that the support of 
the s-equilibrium measure on a torus is, for sufficiently small positive s, likewise a 
proper subset of the torus. 

In this paper we provide sufficient conditions under which the support of the 
equilibrium measure fis,r(A) is a proper subset of the outer boundary of r(A). 
More specifically we show the following. 

• Using rotational symmetry, we demonstrate how to reduce the problem of 
finding the support of the equilibrium measure /Us,r(A) on T{A) for the (sin- 
gular) kernel A;s(x) = l/|x|* to the problem of finding the support of the 
equilibrium measure A^,^ on A for a related kernel ICs which is continu- 
ous when < s < 1 and is singular when s > 1. Lemma 2.2 summarizes 
properties of the kernel /Cg. We further discuss the asymptotics of optimal 
/Cg-energy point configurations on A in both the continuous and singular 
cases. 

• We show that there arc infinite compact sets A for which the support of 
the equilibrium measure on T{A) is all of r(^) for every < s < 1. For 
example, this holds for compact subsets A of a horizontal or a vertical line- 
segment (see Corollary 3.4). 

• We construct sets of revolution r{A) such that the support of the equilib- 
rium measure on T{A) is a proper subset of the outer bomidary of T(A), in 
contrast to the Coulomb case s = 1. We demonstrate this for < s < 1/3. 
(This follows from Theorem 3.7.) An example is the outer boundary of the 
"washer" T{A), where A is the rectangle with lower left corner 1/2 — i/2 
and upper right corner 1 -|-i/2 (cf. Example 3.5). We conjecture that there 
exists for every < s < 1 a compact set A for which supp iJ.s,r(A) is a 
proper subset of the outer boundary ofT{A). 

• We show that for certain sets A the support of the limit distribution on 

sets of revolution T{R + A), for the translate R + A — {R + z | z e A}, 
tends to the full outer boundary of r(i?-|-^) as i? ^ oo. For example, this 
property holds if the outer boundary of A is a compact subset of a circle 
with radius r centered at a > r and < s < 1 (cf. Lemma 5.3). 

• We also show that the support of the equilibrium measure for the loga- 
rithmic case (s = 0) can differ significantly from the case s > 0. For 
example, let ^ be a horizontal line-segment in H+. Then we show that 
supp Ha,T{A) = r(^) for all < s < 1, while it is known that suppy^o,r(^) is 
the circle generated by the "right-most" point of A. (For further discussion, 
sec end of this section.) 

Outline of the paper. In Section 2 we reduce the equilibrium problem to a 
minimal energy problem in the plane with respect to a new kernel Kg for which we 
find an explicit expression. 

Section 3 is devoted to the study of suppAg^A for the kernel /C^. A convexity 
argument (Theorem 3.1) yields that compact subsets A of horizontal or vertical 
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line- segments are examples with supp Xs,a — A for every < s < 1 (Corollary 
3.4). In contrast, we prove the existence of compact sets A for which supp A^, a 
is not all of the outer boundary of A by using the variational inequalities for ICg- 
The essential result here is the 3-point Theorem (Theorem 3.7) which provides a 
sufficient condition for a point on the outer boundary to not belong to the support 
of the equilibrium measure corresponding to /Cg. 

In Section 4 we study the /Cs-equilibrium measure on sets obtained by translating 
a given set A C IHI+ a distance R units to the right. The asymptotic expansion of 
ICs{R + z,R + w), z,w e A, as R becomes large, is given in Lemma 4.1 and it is 
sensitive to the order of the limit processes s — > 0+ and i? — > oo. The relation 
between the energy problem for ICs on A and the energy problem for ICs on the 
translate ii + ^ is discussed. 

In Section 5 we study the kernel that arises as i? — > oo, namely 

We show that any compact subset A of a line-segment [z', z"\ C IHI+ has the property 
that supp = A for every < s < 1, where A^^ defines the equilibrium measure 
for this kernel, and we find an explicit expression for the equilibrium measure A^^ 
on A = [2', z"\. In case that the outer boundary 5 of A is a subset of a circle C we 
get supp A^y^ = S for every < s < 1; see Lemma 5.3. In particular, if 5 = C, the 
equilibrium measure on A for the infinity kernel is simply the normalized arc-length 
measure on C. 

In Section 6 we discuss the discrete Riesz s-energy problem on r(A) C M"^ as 
well as the discrete /C-energy problem on A C 1HI+ for the kernel K. = Kg, K. = ici^\ 
and /C = lci°°K We consider the potential theoretical case < s < dimr(j4) and 
the hypcrsingular case s > dimr(A). In the hypcrsingular case the discrete energy 
problem becomes a weighted energy problem which allows us to use results from 
[4]. We find the limit distribution of minimal /C-energy A/'-point systems, consider 
the separation of such optimal point configurations, and give asymptotics for the 
discrete minimal energy as N 00. Also included are numerical experiments 
showing minimal energy point configurations on Cassinian ovals, line-segments, 
and circles. 

An appendix to the paper provides the computations showing convexity of the 
kernel ICg on the vertical line-segment. 

2. Reduction to the plane, the kernel ICs 

First we fix some notation. The axis of revolution is identified with the y-axis in 
M^. Any vertical cutting plane gives a cross-section of the set of revolution and may 
serve as a reference plane. Selecting a vertical cutting plane we choose one of the 
two closed halfplanes and call it 11+ . It may be identified with the complex right 
half-plane. Then the set of revolution generated by A c M'^:={x+iy \ x >0,y gR} 
is the set 

(2.1) r(A):={R^x|xe A,0<^<27r}, 

where is a rotation by angle (f> about the axis of revolution. The set r{A) is 
obtained by revolving A around the vertical axis. Thus, a single point x + iy € H+, 
a; > 0, becomes a horizontal circle with center on the vertical axis. 
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A Borel measure fi G A^(M ) is rotationally symmetric about the y-axis if 

(2.2) /XR^S) = fXB) 

for all Borel sets B C and for all rotations R<^ about the y-axis. (Here R^B 
denotes the pointwise rotated set {R^x | x e B}.) 

If /i G A^(R^) is rotationally symmetric about the y-axis, then fi can be written 
as a product of two measures, the normalized Lebesgue measure on the half-open 
interval [0, 27r) and a measure jj, on H+, that is 

(2.3) d/t=^dM, = A°r G A1(H+). 

ZTT 

Then the energy of the (compactly supported) measure fi can be expressed as 
Is [A] = / / ks (x - y) d fi{x) d fi{y) 

J JR^xR^ 

J Jm+xu+ 



(2.4) 



where the kernel /Cg {z, w) is given by the integral 

(2.5) ICs{z,w):=^ 



2-K ^ 



|R0^ - w\ 

2.1. The energy problem for /C,. Let < s < 1. Let A C H+ be a compact 

set such that cap, T{A) > 0. Then the uniqueness of the equilibrium measure 
l^s,r(A) oil r(A) and the symmetry of the revolved set T{A) imply that fJ,s,r(A) is 
rotationally symmetric about the y-axis and so dfis.r{A) = [d</'/(27r)] d A^,^, where 
K,A = Ms,r(A) or G M{M+). Furthermore, if G M(M+), then dz>:=[d(/)/(27r)] dz/ 
is rotationally symmetric about the y-axis and so we have 

Js[>^s,a]> inf JaM= inf IsM > Is[iJ-s,a] = Js[K,a]- 
veM{M+) ueMiM.+) 

In the case < s < 1 the equilibrium measure on T{A) is concentrated on the outer 
boundary of r(A) (cf. [12, Ch. II, no. 13]). 

Proposition 2.1. Let < s < dimr(A). Let A C H+ be a compact set with 
cap, r(^) > 0. Then As, a = f^s.r(A) ° T uniquely minimizes Js\y\ over all measures 
v G M{A). Thus, Xs^A is the equilibrium measure on A for the kernel /Cg. It is 
supported on the outer boundary of A. 

The /Cg-energy of a measure was defined in (2.4). The energy V^c^ of A is given 

by 

(2.6) VKSA).= i^i{Js[i^]\y &M{A)}. 
The following relations hold: 

(2.7) VkAA) = Js[\s.a]=IsWt(a)\ = Vsin^)). 
For v G M{A), we define the /Cs-potential by 

(2.8) {z):= ( Ks{z,w) dv{w), ^ G H+. 

J A 
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Let v G M.{T{^A)) be rotationally symmetric with Av ~ [d 0/(27r)] d i/, where v = 
i> o r G M{A). Then the potential Ug is constant on circles r({2;}), z G H+. 
Abusing notation, there holds the following connecting formula 



(2.9) 



U^{z)= [ k,{z-y)Av{y) = ^ ( r K{z - R4,w) d^du{w) 
JriA) Jo 



= / K:s{z,w)du{w) = W^{z), zem+. 

J A 



Prom the properties (1.2) and (1.3) of the equilibrium potential Us^'^'^^^ we infer 
the variational inequalities for Kg for compact sets A in the interior of 13+: 

(2.10) W^'''^>VkS^) everywhere on a, 

(2.11) W^''^<Vk.{A) onsuppA,,A. 

In this case we do no longer need an "approximately everywhere" exceptional set, 
since each point of A generates a circle in with positive capacity. 

2.2. Properties of the kernel /Cg. Let z:=x + \.y, w:=u + w, where x,y,u,v G M. 
Let Wi,:= — w = —u + iv denote the reflection of w in the imaginary axis. 

Lemma 2.2. Let s > 0. The kernel ICg :M+ xB+ in (2.5) has the following 

properties: 

(1) ICs{z,w) is well defined for z ^ w for all s > 0. 

(2) ICs is symmetric: lCs{z,w) = ICs{vj,z). 

(3) ICs is homogeneous: ICs{rz,rw) = r^^K,s{z,w) for all r > 0. 

(4) K.S is continuous at all points (z, w) G H+ x H+ with z ^ w. If < s < 1, 
then Ks is continuous at {w, w) with Re[w] > 0. ICs{z, w) is singular at 
z — w for s > 1 . 

(5) If w is on the imaginary axis and s > 0, then K.s{z, w) = \z~ w\~'^ , z ^ w. 
//Re[w;] > 0, then, for s> 1, the following limit holds: 

(2.12) |^-^r-^/C,(^,t.)^ ^^}!r/y^\ ^ asz^w. 

V7rr(s/2) \w-Wt,\ 

(6) lCs{u + \t,u + w) decreases along vertical lines as \t — v\ grows and lCg{u + 
iy, w + i + iy) decreases along horizontal lines as t > grows^. 

(7) Let < s < 1. For fixed w with Ke[w] > 0, the function ICs{z,w) has 
exactly one global maximum at z = w in BI+. At {w,w) or (w*,w), the 
kernel Kg takes the value 

(2.13) lCs{w,w)=lCs{w,,w) =Xg (s';^) IRewf, 
where 

(2.14) X. fs'; . 2-.^M, = Jll^. 



' 27r; 0Fr(l-s/2) [r(l-s/2)]' 



This follows from differentiating the integral (2.23) with respect to t. 



THE SUPPORT OF THE LIMIT DISTRIBUTION 



7 



(8) The kernel ICs has the following representations in terms of hypergeometric 
functions [2] or in terms of a Legendre function [1] 

(2.15) /C,(0,«;) = |z-«;,r%Fi(«/2.V2;i ' ' 



\z — w] 



I |2 



(2.16) = ^ . .F, ( V2,.A 

[\z - w»\ + \z - w\ ) \ 1 {\z-W:t\ + \z-w\) I 

Observe that the Legendre function is evaluated at values > 1 «/ Re ^ > 
or Re w > 0.^ For s> 1 one can factor out the singularity at z = w, 

(2.18) ^.(.,H = 4^.FJl-Al/2;i_tHlY 

(9) As s ^ Q'^ we recover the logarithmic kernel Kq studied in [11]; 

(2.19) lim ^-("'")-^ =log 2 



0+ S \z — w\ + \z — Wi,] 

(10) As s ^ \~ , lCs{z,w) ^ lCi{z,w), where 

(2.20) K,{zM-=- 1 1 K f | |"-"*|-|^-"| 

IT \z — w*\ + \z — w\ \y\z — Wt,\ + \z — w\ 

and K denotes the complete Elliptic integral of the first kind [1] . 

Remark 2.3. For the special case of the sphere the formula (2.16) reduces to the 
formula (4.14) in Dragnev and SafF [8]. 

The level sets of /Cs(., w), w e H+ fixed, look like Cassinian ovals, cf. Figure 1. 
The asymptotical behavior oiJCs{R + z,R-\-w) as — > oo is given in Lemma 4.1. 



Figure 1. Level sets for /Cs(z, 1), s = 1/2. 



Proof of Lemma 2.2. Let z,w £ EI+ with z = x + iy and w = u + 'w. The relation 

|R0Z — = + V? — 2xucos(f) + {y — w)^ 
gives l/jR^z — = {E ~ F coscj))^'^/'^ for the integrand in (2.5), where we define 

(2.21) E:=x^ + + iy - vf , F:=2xu. 

By (2.21) the kernel Ks{z,w) is symmetric in z,w. The substitution (j) = ip + tt 
yields 

(2.22) Ks{z,w) = ^[ {E + Fcosip)~'''^dil) = - I {E + F cos(j))~'^^ dcj). 

271" J-TT Jo 

Applying the half angle formula and substituting ip ^ (l)/2 we obtain 

(2.23) Ks{z,w) = {E + F)-''^-j^ (^l___sin2^j dV- 
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The integral in (2.23) resembles that of a complete elliptic integral. Indeed, for 
s = 1 this integral is the complete Elliptic integral of the first kind K(fc^) with 
elliptic modulus P = 2F/{E + F). (See for example [1, 17.2.19,17.3.1].) The 
transformation — tt in (2.22) gives 

(2.24) lCs{z,w)=E-'''^- j (^l--cosVj dV- 

The integral in (2.24) is a generalization of Epstein and Hubbells elliptic integral. 

We refer to [16] for a discussion of these elliptic-type integrals. 
A change of variables t = sin^ tjj in (2.23) yields 

(2.25) K,{z, w) = {E + Fr'^ 1 iV2-i (i _ ^f^-^ (^^ _ ^l_t^ ^ t_ 
Recall, that the Gauss hypergeometric series [1, 15.1.1] 

V f! = r(c) ^ Via + k) r{b + k) 



represents the Gauss hypergeometric function 2F1 {^^J^', for all complex z within 

the circle of convergence, the unit circle \z\ = 1. The analytic continuation in the 
^-plane cut along the segment [l,cx3], [1, 15.3.1], 

^^""^ f t^-^il-tf-^-^il-ztydt, Rec>Re6>0, 
Jo 



(2.26) Ks{z, w) = {E + F)-''^ 2F1 (^V2, 1/2. 



T{b)T{c-b) 

can be used to derive a hypergeometric function representation of the kernel JCs {z,w), 

2F 

eT~f, 

Let Wi,:= — w = —u + iv denote the reflection of w in the imaginary axis. Then 

E-F={x-uf + {y-vf = \z-w\\ 

E + F={x + uf + {y-vf = \z-w,\\ 
and we get the relations 

o<-^= '"-"-''-'r"'' = ,<i, z,v:em\ 

- E + F {x + uf + {y-vf~ 

Substitution of (2.27) into (2.26) yields (2.15). 

The hypergeometric function in (2.15) is of the form 2F1 (^'^^i ■ The quadratic 

transformation [1, 15.3.17] yields a more symmetrical representation (2.16). 

In the argument of the hypergeometric function in (2.16) appears the expression 



wJ — \z — w\ \z — w^\^ — \z — 4Re2; Rei 



'\z-w^\ + \z-w\ {\z-w^\ + \z-w\f {\z-w^\ + \z-w\f' 

It satisfies < 1 and equality holds iov z = w ov z = w* only. Therefore we may 
use the series expansion of the hypergeometric function to get 
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If < s < 1, this series converges even for z = w^w*. For z = w, the argument 
of the hypergeometric function in (2.15), (2.16) is 1. Prom [1, 15.1.20] 

r(i-s) r((i-s)/2) 



ICs{w,w) =/Cs («;*,«;) 



IRewl 



/^r(l-s/2) 



iRewl 



m-s/2)Y 

(The first two relations follow from (2.16), the last one from (2.15).) Note, the 
leading coefficient at the right-most is the energy J5(Si;d(/)/(27r)), where is the 
unit circle and d^/(27r) the imiform measure on S^. This shows (2.13) and (2.14). 

Those hypergeometric functions that allow a quadratic transformation are con- 
nected with Legendre functions. From (2.15) and relation [1, 15.4.7] we get (2.17). 

From (2.15) and relation [1, 15.3.3] we get (2.18). From (2.18) follows (2.12). If 
Re[u;] = 0, then w = w*. Hence, by (2.15), ICs{z, w) = \z — w\~^, z w, ioi s > 0. 

The complete Elliptic integral of the first kindK{k'^) [1, 17.3.1] can be represented 
through a hypergeometric function [1, 17.3.9], 

rTr/2 



K(fc" 







^l-fc2 (sini?)^ 



Thus (2.20) follows from (2.16). 

As s — > 0+, the hypergeometric series in (2.28) reduces to 1. Thus it makes sense 



to consider the quotient {ICs{z,w) 
Then 



l)/s. Fix z,w em+ in (2.28). Let z ^ w. 



—JCs{z,w) 
as 



lCs{z,w)log 



+ 



z — w\ + \z 
2 



w\ 



' _d_ 

d^ 



E 

.£=0 



We are only interested in d/Cs(z,u')/ds at s = O"*". }Cs{z,w) becomes one at s = 0. 
The derivative in the right-most term above exists and vanishes. This follows from 



1 



E 



(,s/2) 



and the limit process s 



£\ 



1 



(1 



s/2), 



{i + iy. 



0+. By the ratio test, the infinite series on the right-hand 
side above is absolutely convergent for |^| < 1 (that is z w) and < s < 1. In 
the case z = w one uses (2.13) instead of (2.28). □ 



3. The support of the equilibrium measure for the kernel ICs 

By Proposition 2.1, the equilibrium measure Xs^a on A for ICg is supported on 
the outer boundary 5 of ^. A convexity argument yields sufficient conditions for 
supp Xs.A = S. Recall that a function / : [a, 6] ^ M is strictly convex on [a, b] if 
/(rx + (1 - T)y) < Tf{x) + T)f{y) for all a < x < y < 6 and < t < 1. 

Theorem 3.1. Let < s < 1 and A he a compact set in the interior o/H+. 

(i) 7/7 : [a, 6] IHI+, a < b, is a simple continuous non-closed curve cov- 
ering the outer boundary S of A, that is S C 7*:= {j{t) \ a<t <b}, and 

ICs{"f{-),"f{t)) is a strictly convex function on the intervals [a,t] and [t,b] 
for each fixed t € [a,b], then there is some closed interval I C [a,b] such 
that supp As, A = 7(7) n S. 
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(ii) If J : [0,6] H+ is a simple continuous closed curve, that is 7(0) = 7(6), 
with S' C 7* and extended periodically by j(t) = j{t + b), a/nd /C., (7(-), j(t)) 
is a strictly convex function on the interval [t, t + b] for each fixed t e [0, b], 
then suppAs^A = S. 

Remark 3.2. Note, that S is only required to be a compact subset of 7*. For 
example, S may be a Cantor subset of 7*. 

Proof of Theorem 3.1. Set A = A^,^ and = W^^. We have supp A C 5 C 7*. 
Suppose G is a component of the complement of supp A in 7*. Now observe, that 
by our assumptions, G always corresponds to a subinterval / of one of the sets 
[a, t], \t, b] or \t, t + b] for 'y{t) G supp A. Two cases are possible: (i) Both boundary 
points of G arc in supp A. Then the equilibrium potential W"^ assumes the value 
J'lC, [A] on the boundary of G and, due to strict convexity of 07 on I, is strictly 
less than this value in the open set G. Since > Jjc, [A] on Ad S, no point of 
G is in A. (ii) At least one boundary point of G is not in supp A. This can only 
happen when 7 is a non-closed curve. Without further assumptions the convexity 
property alone is insufficient to show G fl A = 0. From (i) follows the existence of 
some closed interval / C [a, b] such that supp A = 7(/) (1 S. If 7 is a closed curve, 
then 7 = [0,6]. □ 

Remark 3.3. In the proof of Theorem 3.1 we use three main properties: (i) The 
kernel is continuous, (ii) supp A^,^ C S, and (iii) the equilibrium potential satisfies 
a variational principle. These properties also hold for ici°°^ introduced in Section 
5. Therefore, Theorem 3.1 can be applied in case of ici°°\ 

Using Theorem 3.1(i) we next show that any compact subset A of a horizontal 
or vertical line-segment satisfies suppA^,^ = A for every < s < 1. We contrast 
this with the logarithmic case, where it is still true that suppAo,A = A in case of 
a vertical line-segment [11, Cor. 1]. However, in case of a horizontal line-segment 
one has that Ao,a is a unit point charge at the right-most point of A [11, Thm. 1]. 

Corollary 3.4. Suppose A is a compact subset of either (a) the horizontal line- 
segment [a + ic,b + ic], < a < b, or (b) the vertical line-segment [R -\-ic,R-\- id], 
R> 0, c < d. Then supp As,a = A for every < s < 1. 

Proof. For (a) consider the parametrization 7(0;) = x + ic, a < x < b. Prom (2.16), 

/C.(7(^),7(«))=^-^2Fi(V2,V2;^^ 

/C.(7(^),7M)=«-SFi(V2,s/2.g^ 

From (x-"-2")" > 0, n > 0, and (a;^")" > 0, n > 1, we get [ICs{"f{x),jiu))]" > 
for X 7^ u and for every < s < 1. Termwise differentiation is justified by uniform 
convergence for ja; — ttj > 6. By Theorem 3.1, suppA^,^ = 7(7) fl A for some 
/ = [a',b'] C [a, 6]. From the series representations above we observe that the 
kernel ICs{'y{x),j{u)) is a strictly increasing function in x for x < u and it is a 
strictly decreasing function in a; for a; > u. Hence, W^'j"* 07 < ^^[As^a] on [a,b] \ I. 
By variational inequality (2.10), I = [a, 6]. 



(V2)„(V2)„„ 



_ (1 



,2n ^—s — 2n 



n=0 



E 

n=0 



(l)n«! 

(^/2)„(V2)„ 

(l)n«! 



, X > u, 



X < u. 
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For (b) consider the parametrization ^{y) = R + iy, c < y < d. From (2.15), 

(3.1) lC.{7{y)Mv))=[^R' + {y-vfY''" '^'{'^\'^' uR4t-vf )- 

A direct calculation (assisted by Mathematica, see Appendix A for more details) 
shows that d'^[K:s{j{y),-f{v))]/ dy^ > for every < s < 1. By Theorem 3.1, 
supp As. ^ = 7(/) n A for some / = [c',d'] C [c,d]. From the representation above 
we observe that the kernel ICs(s/{y),j{v)) is a strictly decreasing function in y for 
growing |j/ — w]. Proceeding as in part (a) we get I = [c,d]. □ 

In contrast to the horizontal or the vertical line-segment we will show that there 

are compact sets A in the interior of IHI+ for which, in fact, the support of the 
equilibrium measure on A for JCg is a proper subset of the outer boundary of A. 

Example 3.5. Let A be the rectangle with lower left corner 1/2 — i/2 and upper 
right corner 1 + i/2. Using Theorem 3.7(c) below with a; = 1/2 and z' = 1 + i/2, 
it follows that 1/2 ^ suppA^ ^ for < s < 1/3. Alternatively, if A is the left- 
half circle with radius 1/2 centered at 1, it again follows from Theorem 3.7(c) that 
1/2 ^ suppAg^A for < s < 1/3. In contrast, as A is moved to the right R units 
and ii —> DO, we get supp A^^^ = A; see Lemma 5.3. 

To prove Theorem 3.7 we use a special case of the following observation. 

Lemma 3.6. Let < s < 1. Suppose A is a compact set in the interior o/H+. 
Let A denote the unique equilibrium measure on A for ICg ■ If 

(3.2) ]Cs{z-,-)> I K.s{-,w')Av{'w') everywhere on snpp X 

Jb 

for some subset B c A and some probability measure v e M.{E), then z ^ supp A. 
Proof. Using (3.2) and the variational inequality (2.10), we get 

W^^(0) = j Ks{z,w)dX{w)> j j Ks{w,w')dv{w') A\{w) 

= I W^^(w;')dKw'')>^;cJA] [ Av{w')=JkA>^]. 
Jb Jb 

But W^^(2;) > jT/c^W implies, by the variational inequahty (2.11), that z ^ supp A. 

□ 

Let z = X > and set B = {z',z'}, z' in the interior of IHI+, Im[z'] ^ 0, and 
place the charge 1/2 at each point in B. Then (3.2) is equivalent to the property 

(3.3) Kls{z,-) > K,*{-,z') everywhere on supp A, 
where /C* denotes the kernel 

(3.4) K* {z, w):= [K, {z,w) + {z, w)] /2. 

Theorem 3.7 (3-point Theorem). Let < s < 1. Let x > and z' he in the 
interior of M.'^ . Let A be a compact subset of {w &'E.'^\lCs{x,w) > ICs{x, z')} in the 
interior 0/H+ with x, z', z' e A. 

(a) If As:=Kls{x,z') — K*{z' ,z') > 0, then x ^ suppAs,^. 
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(b) // z' = 1 + i7, 7 > 0, and condition 

(3.5) 4 (7 + yr+T^) > (^7(1+^7+7^ + ^(l-a;)'+72^ 

is satisfied, then A,, > (and hence, by (a), x ^ suppA^^/ij for s > 
sufficiently small. 

(c) If X = 1/2 and z' — l + i/2, then A^ > (and hence, by (a), x ^ suppAg^A^ 
for alio < s < 1/3. (The graph of Ag is shown in Figure 2.) 

Proof of Theorem 3.7. We show first (a). The function /Cs(z, •) has a unique max- 
imum at z in H"*" (Lemma 2.2(7)). So 

(3.6) JCsix, w) - Kliw, z') > Ks{x, w) - K.*{z', z') > lCs{x, z') - Kl{z', z'). 

The first inequahty holds in IHI+. The last one holds on {w G M.'^\lCg{x,w) > 
ICs{x,z')}. Now, let A be a compact subset of {w G M.'^\)Cs{x,w) > )Cs{x,z')} in 
the interior of EI+ with ' e A. Then 

W:{x) > [W^{z') + W^i^)] /2 + /C,(x, z') - Kiz', z'), V e M{A). 

This follows from (3.6) and ICs{w,z) = K,s{z,w). If the difference As:=/C,s(.x, z') — 

lC*{z', z') is positive, the variational inequahty (2.10) imphcs Wg"'* {x) > J/c^ [^s,a]- 
Therefore, x (/ suppA^./i, by variational inequality (2.11). This shows (a). 

Set z' = /3 + 17 with /?, 7 > 0. Since ICs{pz, pw) = p~^lCs{z,w), p > 0, we may 
fix one of the variables x, /3, or 7. Let /? = 1. From (2.15), (3.4), and Lemma 2.2(7) 
we get 



(3.7) 



{l + xf+j' 



^ / s/2. 1/2. 4.x 

1 '(l+.x)V7^ 



2 



- '-2-^ (1 + 7^) 2F, fV2, 1/2. _J_) _ i2-«ifci)Z^. 
2 ^^^^ ) ^^n^ ^ '1+72 2 V^r(l-s/2) 

We approximate Ag by its series expansion at s = 0. Prom Lemma 2.2(9) and (3.4) 
As 

lim — = /Co(x, 1 + 17) - /Co(l + i7, 1 + 17) 

s->0+ S 



4(7+v/T+^) 
(^^{l + xf+'y^ + ^{l-xf+^^^ 



\ log ' ■ ' . - > 



which implies (b). 

We show that A^ (as a function in s) is strictly concave on (0, 1) if a; = 7 = 1/2. 
Using (3.4) and integral representation (2.25) we get 



A,:=/C,(x,l + i7)-/C:(l + i7,l + i7) = ^ / t-"''{l-t) ''^g{s,t)dt, 



where 



g{s, t):= [(1 + xf + 7^ - 4xt\ -^[4(1+ 7') - 4i] - 1 [4 - 4*]"^/' 
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Negativity of {d/ds)^g{s,t) for all < i < 1 implies that Ag is strictly concave. 
Let a; = 7 = 1/2. Prom {d/dsfr-^l'^ = (l/4)F(s,r), F{s,r):=r-^l'^{\ogrf , we get 

4 g{s, t) = F{s, 5/2 - 2t) - {l/2)F{s, 5 - 4t) - {l/2)F{s, 4 - At). 



Negativity of the right-hand side above is equivalent with 

'5/2 -2ty/' [log(5-4t)]2 /5/2-2ty/' [\og{A-4t)f 



2I+S/2 ^ 



5/2 -2t J [log (5/2 -2f)]' V 2-2t / [log (5/2 - 2^)]' ' 



For growing s the left-hand side of the last relation is increasing while the right-hand 
side above is decreasing. Thus 

(3.8) 2^^V. < J^^il^^ + [^o^i^-^t)r^ 
[log(5/2-2t)f [log(5/2-2t)f 

is a sufficient (but not necessary) condition for {d / ds)^ g{s,t) < for all < i < 1. 
Elementary calculus shows that hi{t) > 3 on [0,3/4) and h2{t) > 4 on [3/4,1). 
Hence the right side of (3.8) is > 3 > 2^^"^ on [0,1). Consequently, is strictly 
concave for < s < 1 (cf. Figure 2). 

Since has a zero at s = with lims^g+ A^/s = (1/2) log(\/5 — 1) > 0, 
limjj^i- As = —00, and A^ is strictly concave on (0,1), the difference Ag has 
exactly one other zero in the interval [0, 1) denoted by si. By our reasoning Ag > 
if and only if < s < si. A numerical solver gives si « 0.341107. . .. Numerical 
computation shows that A1/3 « 0.0011 > 0. This can be rigorously justified by 
assistance of Mathematica and use of exact arithmetic. □ 

By Theorem 3.7(a), the positivity of Ag implies x ^ suppAg^^- By (3.7), Ag 
depends on three parameters x, 7, and s. Sec Figure 2 for a plot of the level surface 
As = 0. This 0-level surface is the boundary of the set of admissible configurations 
{x, 1/7, s) using a three point scheme = a;, 2;' = 1 -|- 17, and z'. Prom Figure 2 we 
get numerical evidence that the maximum s possible for a three point approach is 
about 0.38. 



Figure 2. 0-level set of A^ in (3.7) cut off at I/7 = 4 and Ag for 
X = 1/2, 7 = 1/2. 



4. Kernel ICs in the limit R ^ 00 
We want to study the behavior of ICs{R + z, R + w) as R becomes large. 
Lemma 4.1. Let < s < 2, s 7^ 1, and z,iu & H+. Then 



V 27r y 1-s V7rr(l + s/2) 2R 

w/iere Xs(S^; d0/(27r)) is given in (2.14). 
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Remark 4.2. In the case 1 < s < 2 the second term in (4.1) becomes the dominant 
term. In the special case s = 1 the following expansion can be shown: 

(4.2) 



2R}Ci{R + z,R + w) = + -logR 



log \z — wl 

n 



R 



00. 



Proof of Lemma 4-1- Let < s < 2, s ^ 1. Using [1, 15.3.6], we get a representa- 
tion of (2.15), 

)Cs{R + z,R + w) 

|2 \ 



^ r(l - s/2) 

2 r((l + s)/2) \z-w\^ 



(l + «)/2' \2R + z 



Fil 



1- .s/2, 1/2 \z-w\- 



1-s ^/i^(s/2) \2R + z-w4' + - s)/2' \2R + z - w.f J' 

with convergent series expansions of both hypcrgcomctric functions. The first one 
is of the form 1 + 0{sR^^), the second one is of the form 1 + 0{R~^). Since 



1 + 



2R 

we get 

K:s{R + z,R + w) 

r((i-,s)/2) 



s Re[z — 
~ 2 W 



+ 



R 



oo, 



= 2" 



R- 



r(i - s/2) 
1 r((i + s)/2) 



s Rc[z — 



\z — w 



2 R 

l-s ^-1 



l + 0{ 



(i?2 



l + O 



l + O 



1 

R^ 



l-s ^/^^T{s/2) 
We reorder the terms with respect to powers of R and obtain (4.1). 

It is convenient to define the following kernels 

(4.3) ICi'^^z,iv):=2R[lCs{R + z,R + w)-Is{^^;<i(j)/{2n))R-'], < s < 1, 

(4.4) ICi^\z,w):=2RICs{R + z,R + w), s > 1, 
and 

2 r((n-.)/2) r((.-i)/2) 



□ 



(4.5) )Ci^\z,w):= . 

1 — S ^/TT 1 (s/2) 

Then, by (4.1), 

(4.6) lim K:i^\z,w) = ICi°°^( 

R — *oo 

and, from (2.18) and [1, 15.1.20], it follows 
(4.7) 



/tF r(s/2) 

< s < 1, 



lim Iz-wl'-^JCf '^z^w) = \z~w\'-^ ICi°°\z,w), 

R^oc 



S > 1, 



where in both cases the convergence is uniform on compact subsets of H+ x H+ . If 
s < dimr(A), we let Jj^(R)\y\ and J^(^-)\y\ denote the associated energies of the 

compactly supported measure v e A^(BI+). From the definition of the kernel /ci^^ 
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we see that the equihbrium measure A^^ on the compact set A C 11+ for the kernel 
/Cs^^ is equal to the equilibrium measure Xs,r+a on R + A for the kernel JCg in the 
following sense: Af^^(B) = XA+niR + B) for a Borel set BcB+. 

Remark 4.3. The asymptotics (4.1) holds uniformly in < s < sq < 1. So 

lim /C(^) {z, w)/s = l6^^ {z, w) + {l/R) , i? ^ oo. 

The expression z—Wi,] — \z—w\ is the oo-kernel for the logarithmic 

case introduced in [11]. However, reversing the order of limit processes, we get 

lim /C(«)(^,«;)/s = /C(~)(^,w;)/s. 

it— »CX) 

Now, in the limit s — > 0+, the right-hand side above tends to —\z — w\. 

5. The energy problem for the kernel /ci°°^ 
5.1. The case < s < 1. The kernel 

falls into a class of kernels studied by Bjorck [3]. From his results wc infer that to 
every compact set A C H+ and every < s < 1 there exists a unique equilibrium 
measure A^^ supported on the outer boundary of A. ( "Outer boundary" is justified 
by the strict superharmonicity of the infinity kernel everywhere in C.) Let W;^, 

denote the potential for a measure /x e M.{A) and for the kernel /Cs°°^ 

w;; {z):= f /c(°°) {z, w) d ^i{w\ ^ e e+. 

'^o J A 

Then W^(^, is continuous on H+ and from results in [3] there follows that W^°'^, > 

i/^Ccx,) [A^^] on A and equality holds on suppA^^. We note, that A^^ converges 
weak-star to A^^ as i? oo. This follows from the weak-star compactness of 
A4{A), relation (4.6), and the uniqueness of the equilibrium measure A^^. 

Suppose the curve 7 : [a, b] —> 11+ covers the outer boundary S of A. Set 
= 1 7(0 — w\. Assuming 7 is twice differentiable at t we have 

(5.1) ^JCi-\,it),.) = 2^^|±^ [. ir'S - r.r: 

Then for fixed u>, wc have that ici°°\j{t),w) is strictly convex on any interval 
where s{r'^,)'^ > 'r„,r'',. A sufficient condition would be r'^ < 0. 

In the following we give examples of compact sets A C IHI+ such that the support 
of the equilibrium measure on A is given by the outer boundary of A. 

Lemma 5.1. Let A be a compact subset of the line-segment [z' , z"\ in the interior 
o/H+, z" -z' = 2re''^, r > 0, < < tt. Then suppA^^ = Afor aU^<s< 1. 
In particular, if A = [z',z"], then 

(5.2) dX-A^o) = ^^ + ;)ff r^- (r^ - T^)^/^"^ dT, 
where w = {z' + z")/2 + Te''^, |T| < r. 
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Proof. W.l.o.g. consider the parametrization 7(t) = te"^, \t\ < r. Then 

^/C(-)(,(0, = 2'-^^ \t - > 0, < . < 1. 

By Theorem 3.1 there exists an interval I C [—1, 1] such that supp = 7(7) fl A. 

Since the kernel lci°°\'j{t),j{T)) decreases as \t ~ T\ grows, there follows that the 
equilibrium potential is strictly less than J^(ac) [A^^] on (— ooe"^, ooe"^) \7(/). But 
the equilibrium potential is > J^(,y=) [A^^] on A. So, supp A^^ = A. Relations (5.2) 
follows from the constancy of the integral 

\t - T|^-" (r^ - T^)"/^-^ dT = r(s/2) r(l - s/2) 



/ 

J — g 



and the fact that the /ci°"'' -potential for this measure (5.2) is strictly decreasing 
away from the line-segment. We used the auxiliary result Lemma 5.2. □ 

Lemma 5.2 ([14, Hilfssatz I]). Let -\ < a <\, ai- 0. Then for -l<y< 1: 



2N,-(l+a)/2 1^ _ la J ^fl-a\_^fl + a 



(l-a;2)-™'^|x-2/rdx = r 



cos(7ra/2) ' 



Lemma 5.3. Let the outer boundary S of the compact set A be a subset of a circle 
C centered at a > with radius < r < a. Then supp A^^ = S for every < s < 1. 
In particular, if S = C, then A^^ is given by the normalized arc-length measure on 
C and supp A^^ = C for all < s <1. 

The result supp A^^ = 5 for < s < 1 differs considerably from the logarithmic 
case. By [11, Thm. 4], one has suppA[f^ = {a + re'"^ | \(p\ < 6} for some 9 e [0,7r/3]. 



Proof. W.l.o.g. consider the parametrization 7(1/)) = re"^, < 4> < 27r. Then 

22-- r((l + s)/2) 



/Ci°°)(7(<^),7(</'')) = -Y- 



sm ■ 



l-s 

l-s 



r 



s 0Fr(s/2) 
By direct calculation (assisted by Mathematica) 

^^(-)(7(0) 7(0')) = ,-s-, m + s)/2)l + s-il-s)oos{c^-<l.') 
^^,l^s (7(0j,7(0)) V^T{s/2) |sin[(<^-<^0/2]r+^ 

Since 7 is a simple closed continuous curve and j{(t>) = j{(t> + 27r), by Theorem 3.1, 
supp A^^ = 5*. In the case S — C, rotational symmetry gives d A^^ = (27r). □ 

5.2. The case s > 1. The kernel 



1 



V^r{s/2) \z- 



-1 ' 



s > 1, 



is (up to a multiplicative constant) the Riesz-(s — l)-kernel in the plane which 
can be identified with C. If 1 < s < l + dimA, then classical potential theory yields 
that there exists a unique equilibrium measure A^^ on A with supp A^^ D A, where 
A denotes the set of all points of A each neighborhood of which intersects A in a set 
of positive Riesz (s — l)-capacity. Examples of sets A with A = A are line-segments, 
circles, or more generally, any Jordan curve; discs, "washers". 

Lemma 5.4. Let A be a compact subset ofC with dim A > and s a real number 
with 1 < s < 1 -|- dim^. Then Xf^ converges weak-star to A^^ as 00. 
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Proof. To simplify notation we use the abbreviations ICr = ici^\ /Coo = K^s°°\ 
Xr = Af^, and X^o = A^^. Prom the definition of )Cr given in (4.4) and the 
formula (2.18) it follows that 

Kr{z,w) = VtR{z,w)Koo{z,w), {z,w) e C X C, 

where 

We remark that Q.r converges uniformly to 1 on compact subsets of C x C as 

R oo. 

Since M.{A) is weak-star-compact, there exists a weak-star cluster point A* of 
A/j as -R — > oo. We will show that Jk^ [A*] < Jk^ [Aoo] from which the Lemma 
will immediately follow. Let i?^, A: > 1, be a sequence of numbers such that 
lim/j^oo-Rfc = oo and Xr^ A A* € M{A) as k ^ oo. Thus riflj,(Aijj, x Xr,^) A 
A* X A* as fc — > oo and we have (see [12, Lemma 0.1]) 

Jjc^ [A*] < liminf JicH, [Xr,] < liminf J)c^, [Aoo], 

fc— *oo k—^oo ^ 

where the second inequality follows since Xr^ minimizes Jkr^ ■ Finally, since ^r 
converges uniformly to 1 on A x ^ as i? ^ oo, we have liminf;;_>oo JiCr^ [Aoo] = 

Jk^[Xoo\ which shows that Jk.^[X*] < J7k;„„[Aoo]- Since Jk:^ has a unique mini- 
mizer, it follows that Aoo is the only weak-star cluster point of A^ as i? — !■ oo. □ 

In the hyper-singular case s > 1 -|- dim A = dimr(y4.), both energy integrals 
Jk^(r)[i^] and J7'^(oo)[i'] are infinite for every v G A4{A). In Section 6.2 and 6.3 

we consider the limit distribution of minimal /C^-energy and /ci°°^ -energy iV-point 
systems as N ^ oo for the hyper-singular case and for sufiiciently "nice" sets A 
(namely d-rectifiable sets). 

6. Discrete Minimum Energy problems on A c IE1I+ 

In this section we discuss the discrete Riesz s-energy problem on T{A) C M.^ as 
well as the discrete /C-energy problem on A c IHI+ for the kernel /C = /Cg, /C = ]Ci^\ 
and K, = ]Ci°°\ The N -point Riesz s-energy ofT{A) is defined as 

£s{T{A),N):=mmE,{XN), 

where the minimum is taken over all iV-point configurations C r(A) and 
Es{Xn) is defined as in (1.1). We let Xj^ = X'^ ^ denote an A^-point configu- 
ration in T{A) attaining this minimum. 

Similarly, for an A''-point configuration Zn = {zi, . . . , Zn} C A, let 

EK{ZN):=J2lCizj,Zk) 

and let the N -point IC-energy of A be defined as 

£K{A,N):=mmEK{ZN), 

over all iV-point configurations Zjv C A. This minimum is attained at a minimal JC- 
energy iV-point system = {zl, . . . , z^}, that is £k:{A, N) = Ejc{Z^). Finally, 

let X{Z'f^):={l/N)J2k=i^zt>' ^® interested in the weak-star convergence of 
X{Z^) and in the asymptotic growth of £ic{A, N) as N ^ oo. 
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6.1. The potential theory case. For < s < dimr(>l), there is a unique equi- 
librium measure \k,a minimizing the /C-energy 



J^[X] = j ]C{z,w)d\{z)d\{w) 



over measures A G M{A). (See Proposition 2.1 for the case JC = JCg and Bjorck [3] 

for the case K, — ici°°\) 

Proposition 6.1. Suppose A is an infinite compact subset in the interior o/]HI+. 
Let K = Ks,K = ]Ci^\ or K = /ci°°' and 0< s < dimr(^). For N >2, let Z*^ 

be a minimal IC-energy configuration of N points {z*, . . . , z^} C A. Then X{Z^) 
converges weak-star to the equilibrium measure Xk.,a on A as N —^ oo and 

(6.1) lnnJj^=JA,,]. 

Proof. The proof follows using standard arguments as in [14], [12, pp. 160-162] and 
[10]. The essential ingredients of the proof arc the boundedness of the sequence 
(6.1) and existence and uniqueness of the equilibrium measure X/c^a- O 

In Figure 3 we show minimal /Cs-energy configurations for A'' = 32 points re- 
stricted to a Cassinian oval for various values of s with < s < 1. A somewhat 
surprising result of these numerical experiments is that for fixed A'' and for s close 
to 1~ a rather large part of the Cassinian oval is free of points. (Note, that in the 
case s = 1, the support of the equilibrium measure Xs^a is A.) In Figure 4 we show 
minimal JCf -energy configurations for A'' = 32 points restricted to a Cassinian oval 
for various values of s with < s < 1. 



Figure 3. Minimum /C^-energy configurations {N = 32 points) 
for s = 0, 0.5, 0.75, and 1. 



Figure 4. Minimum /C^-energy configurations {N = 32 points) 
for 8 = 0, 0.5, 0.75, and 1. 



Numerical experiments for a circle C centered at a > with radius r {0 <r < a) 

suggest that for a fixed s, say s = 1/4, the equilibrium measure A, ^7 is concentrated 
on a proper subset of C. (For what we can prove, see Example 3.5.) However, the 
equilibrium measure A^,^ associated with the translate R + C converges weak-star 
to A^^ as i? ^ 00 and wc can show (sec Lemma 5.3) that A^^^ is the uniform 
measure on the circle C for all < s < 1. This phenomenon that the support of 
Af seems to spread out as i? — > 00 is illustrated by considering discrete minimal 
/Cs-encrgy points on the translate R + C for varying values of R. In Figure 5 we 
show minimal /Cs -energy configurations for iV = 40 points restricted to translates 
i? + C of the unit circle C centered at a = 1, where R = lO*'/^ (fc = 0, 1, . . . , 5). 
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Figure 5. Minimum /Cg-energy configurations (TV = 40 points) 
on translates i? + C of the unit circle C centered at 1 for R = 10*^/^ 
(fc = 0,l,...,5). 



6.2. The hypersingular case for d-rectifiable sets. Suppose s > dimT{A). In 
this subsection we require that A he a. c?-rectifiable subset of the interior of 1HI+ 
{d = 1,2). Recall that a set X C 11'^ is d-rectifiable, d < p, \i it is the image of 
a bounded set B in M"^ with respect to a Lipschitz mapping, that is a mapping 
: B — > that satisfies for some positive constant c 

(6.2) \(j){x)-(l){y)\<c\x-y\ for all a;, y e S. 

Note that every compact subset of is 2-rectifiable. Also note that if A is a 
rf-rectifiable set in H+, then T{A) is a (d + l)-roctifiablo set in R^, d = 1,2. In 
order to avoid complications, we require that A is in the interior of H+. In this 
case dimr(A) = 1 + dim A. 

Using the properties of ICg as given in Lemma 2.2 it follows that ICs{z,w) = 
fl{z,'w)\z — where : A x A ^ R is continuous and positive. In the 

terminology of [4] then O is a GPD weight function on A (see [4] for the general 
definition of CPD weight function). Also, note that 

^ ' ' V^r(s/2) \w-w^ 

If A is a compact set in R'' and is a CPD-weight function on A x A, then for 
s > d one can define the weighted Hausdorff measure nf on Borel sets B <zM+ 
by 

(6.3) n'f{B)~ [ [n{w,w)]-''^'dndiw). 

JBnA 

Then the following result is a corollary of Theorem 2 in [4] . 

Proposition 6.2. Let d = 1 or 2 and suppose A is a compact d-rectifiable set con- 
tained in the interior 0/H+ with positive d- dimensional Hausdorff measure HdiA). 
Let s > dimr(A). For N > 2, let be a minimal ICg-energy configuration of 
N points {zl'^ , . . . , z^i^ } C A. Then the sequence Z'^, N > 2, is asymptotically 
uniformly distributed with respect to ; that is, 

, N o/s-i,n 

(6.4) 4 y S.^i.- A , as 00. 

Moreover, the minimal N -point Kg-energy satisfies 

(6.5) hm - ^-^"^ 



where Cs~i.d is a positive constant which does not depend on A and N . 

Remark 6.3. The constant Cs-i,d is exactly the same constant which appears in the 
analogue of (6.5) for the non-weighted case, that is for Q.{z, w) = 1 for all z,w & A. 
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It can be represented using the Ricsz s-cncrgy for the unit cube in Mf^ via 

It was shown in [13] that Cs^i = 2^(s), where ({s) is the classical Riemann zeta 
function. However, for other values of d, the constant Cs,d is as yet unknown. 

In the boundary case .s = dimr(^) and for a 1-rcctifiablc set A an additional 
regularity condition is needed to prove a result analogous to Proposition 6.2. The 
following result is a corollary of Theorem 3 in [4] . 

Proposition 6.4. Let d = 1 or 2 and suppose A is a compact d-rectifiable set 

contained in the interior of H"*" with positive d- dimensional Hausdorff measure 
Tid{A). If d = 1, we further require that A is a subset of a -curve. Let 
s = 1 + d. For N > 2, let be a minimal )Cd+i-energy configuration of N 
points {zf~^^'^ , . . . , z'}^^''^} C A. Then the sequence Z'^, N >2, is asymptotically 
uniformly distributed with respect to H^' ; that is, 

1 ^ j^d,n 
— ySd+i,N^ — , as N ^ oo. 

Moreover, the minimal N -point Kg-energy satisfies 

£k,UA,N) 



(6.7) lim 



where (3d = 7r''/^/r(l + d/2) is the volume of the unit ball in W^. 

Remark 6.5. It is a consequence of Theorem 2 and 3 in [4] that minimal Riesz 

s-energy point configurations X'^ C r(^) are asymptotically uniformly distributed 
with respect to Hd+i restricted to T{A) in the hypersingular case s > d + 1. For 
z e 1HI+, let 6z denote the rotationally symmetric probability measure supported on 
r({z}). When s < d + 1 we have that both (1/iV) Ea;GJf* i'^/^)J2zez* 
converge weak-star to Hs,r{A)- However, for s > d + 1, the discrete probability 
measure (1/-/V) J2xex^ converge weak-star to Hd+i (normalized and restricted 
to T{A)), while Proposition 6.2 implies that (l/N) J2zez'^ converges to a measure 
that depends on s. In the boundary case s = d + 1, the latter limit distributions 
are equal (cf. Proposition 6.4). 

In the following we consider two examples: a line-segment in general position 
and a circle centered on the real axis. 

Example 6.6. Let A be the line segment with parametrization j{t) = R + fe"^, 
\t\ < 1, where R > cos(p and < </> < tt fixed. Then A is a 1-rectifiable set and 
the weighted Hausdorff measure can be explicitly calculated. Indeed, since 

dHi(t) = di, we get for s > 2 

ffi8l ^r^'" (,y- (R + tcoscbf^'-'^dt 

^'•'^ Hr^-"(^) " ^ - 1 {R + cos^f'-'^ - (i?-cos0)-/^-^) ' "- 

s cos 6 
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Note, in the case of the vertical hne-segment (that is (j) = 7r/2), the last expression 
reduces to 

j-t I s — 1,0 

(6.9) d ^pf^ (t) = -dt, |t|<l;s>2. 

In Figure 6 we show minimal /C2 and /C4-energy configurations for = 40 points 
restricted to the line-segment with i? = 3/2 and ^ = 7r/4. 



Figure 6. Minimum /Cg-energy configurations (TV = 40 points) 
on a line-segment for s = 2 (left) and s = 4 (right). 



Example 6.7. Let A be the unit circle centered at i? > 1. Then A is a 1-rectifiable 
set and the weighted Hausdorff n 
d Wi (0) = d (/>, one has for s > 2 



set and the weighted Hausdorff measure ^'^ can be explicitly calculated. Since 



(6.10) d^Ho (<^) = \ ' -7r<</.<7r. 



In Figure 6.2 we show minimal /C2 and /C4-energy configurations for A'' = 40 points 
restricted to the unit circle centered at i? = 3/2. 

caption Minimum /Cs-encrgy configurations (TV = 40 points) on the unit circle 
centered at 3/2 for s = 2 (left) and s = 4 (right). 

Remark 6.8. Results similar to Proposition 6.2 and Proposition 6.4 also hold for 
the kernel In this case the diagonal of the CPD weight function becomes 

^ r((g-i)/2) w^w. 

U {w,w):— — ;= — - — --— 1 — 

In particular, wc have the limit 

lim Q''{z,w) = ^^)fr,^y,f , s > dimr(A), 

where the convergence is uniform on compact subsets of H+ x EI+. Consequently, 

o/s-i,n^ -Tj I 
(6-11) T7^i5 -^^TT^, as i?^ 00 and s > dimr(A). 

Here 'Hd\A/'Hd{A) is the limit distribution of minimal A^^-energy iV-point config- 
urations as ^ 00 (see next subsection). 

Of interest is the question of how well minimal /C-energy points are separated, 
that is, we are asking for a lower bound for the separation radius 

(6.12) 5{Z*j^):=m:m.{\z-w\ \ z,w & Zl,,z ^ w} 

of optimal /C-cncrgy A^-point systems valid for N > 2. In fact, such an estimate 
can be obtained on sets of arbitrary Hausdorff dimension a. The following result 
is a corollary of Theorem 4 in [4] . 
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Proposition 6.9. Let < q < 2. Suppose A is a eompaet subset in the interior 
o/e+ with l-Lc{A) > 0. LetJC = JCs or JC = /Cf tinth R>0. Then for every s>a 
there is a constant Cs = Ca{A,Q,a) > 0, where is the CPD-weight function asso- 
ciated with K, such that any K-energy minimizing configuration Z'^ on A satisfies 
the inequality 



(6.13) 6{Z*^)>I ' , ' N>2. 

^ ' ^ - ^ c„(iVlogAr)-i/" s = a, 



6.3. The hypersingular case for the kernel JC^. Suppose s > 1 + dimA. The 

kernel /C^ can be written as ]Cf{z,w) = VL°°{z,w)\z — w\^~'^, where the CPD 
weight function il°°{z,w) = r((s — l)/2)/[-0rr(s/2)] is a positive constant. Thus, 
we can apply the theory developed in [4] to obtain 

Proposition 6.10. Let d = 1 or d = 2. Suppose A is a compact d-rectifiable 
set contained in the interior 0/H+ with positive d-dimensional Hausdorff measure 
TidiA). Let s > 1 + dim A. For N > 2, let he a minimal KLf -energy con- 
figuration of N points {z^^ , . . . , z^i^ } C A. Then the sequence Z'^, N > 2, is 
asymptotically uniformly distributed with respect to Hdi that is, 

Moreover, the minimal N -point Ks-energy satisfies 



(6.14) hm 



jV^oo A/'l+(s-l)/<i [Wrf(yl)]^*"^^/'^' 

where C"^-^ d positive constant which does not depend on A and N . In fact, 
C's^i d = Cg-i d^{{s — l)/2)/[V7r r(s/2)] and Cg-id is the same constant as in 
(6.6)! 

Remark 6.11. The first part of Proposition 6.10 holds for the boundary case s = 
l-\-d, d = dim A, as well. (In the case d = 1 it is also required that A is contained 
in a C -curve.) The mininial iV-point Kis~ 

energy satisfies 



(6.15) lim 



£k,+M,N) _ (3d 



N^oo N'^logN nj{A)' 

where /3d = n^/'^/ T{1 + d/2) is the volume of the unit ball in R^. 

We remark that a separation result like Proposition 6.9 can also be stated for 
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Appendix A. Convexity of the /Cs-kernel on vertical line-segments 
The parameter u is fixed. The second derivative of (3.1) with respect to y is 



where A denotes the difference {y — v). Applying to each hypergeometric function 
the Hnear transformation [1, 15.3.3] and simphfying we get 
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